Compound nuclear decay and the liquid to vapor phase transition: a physical picture 
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Analyses of multifragmentation in terms of the Fisher droplet model (FDM) and the associated 
construction of a nuclear phase diagram bring forth the problem of the actual existence of the 
nuclear vapor phase and the meaning of its associated pressure. We present here a physical picture 
of fragment production from excited nuclei that solves this problem and establishes the relationship 
between the FDM and the standard compound nucleus decay rate for rare particles emitted in 
first-chance decay. The compound thermal emission picture is formally equivalent to a FDM-like 
equilibrium description and avoids the problem of the vapor while also explaining the observation 
of Boltzmann-like distribution of emission times. In this picture a simple Fermi gas thermometric 
relation is naturally justified and verified in the fragment yields and time scales. Low energy 
compound nucleus fragment yields scale according to the FDM and lead to an estimate of the 
infinite symmetric nuclear matter critical temperature between 18 and 27 MeV depending on the 
choice of the surface energy coefficient of nuclear matter. 

PACS numbers: 24.10.Pa,25.70.Pq 
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After decades of theoretical and experimental studies, 
recent papers have presented what can be considered a 
quantitative, credible liquid-vapor phase diagram con- 
taining the coexistence line up to the critical temperature 
for small nuclear systems |l|, |2( . This diagram was not 
obtained by "traditional" methods of caloric curves 0, 0] 
or anomalous heat capacities 0. _]. Rather, it was gen- 
erated from the fitting of the charge distributions in mul- 
tifragmentation by means of a Coulomb corrected Fisher 
droplet model (FDM) 0,00,01 which gives the cluster 
composition of a vapor: 



(rp s A -r MM c eA a 
n A(T) = q A exp I — — 



(1) 



where: go is the normalization, A is the cluster mass num- 
ber, r is a topological critical exponent, cq is the surface 
energy coefficient, T is the temperature, e = (T c — T)/T c 
is the relative deviation from the critical temperature T c , 
and <7 is the surface to volume exponent. The quantity 
cqeA" is the surface free energy cost of forming an A- 
sized cluster. A/i is the difference of chemical potentials 
between the liquid and the vapor. Finite size effects 0], 
Coulomb energies [lfl | and other nuclear energy terms 
(e.g. due to isospin and angular momentum) can be ab- 
sorbed into the chemical potential. In references 0] 
this was done for the Coulomb energy by writing A/i 
as the Coulomb energy of a residue plus fragment "sad- 
dle" configuration compared with the initial state. For 
A/i = the liquid and the vapor are in equilibrium and 
Eq. (JIJ can be taken to be the equivalent of the coexis- 
tence line 00. 

In fact, one can immediately obtain from Eq. Q the 
more conventional p, T and p, T phase diagrams by re- 
calling that in the FDM the clusterization is assumed to 
exhaust all the non-idealities of the gas. It then becomes 
an ideal gas of clusters. Consequently, the total pressure 



and density can be immediately calculated 
p(T)=TY / n A (T), 

A 

p{T)=Y J An A {T), 
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and 



as well as the corresponding scaled quantities p/p c 

P/Pc _____ 

Tests on the lattice gas (Ising) model liilllj,ll2| demon- 
strate a beautiful agreement between the vapor cluster 
concentrations and Eq. JJ|, and analysis of many multi- 
fragmentation reactions JlJ, |2| show equally good agree- 
ment, leading to a characterization of the liquid-vapor 
phase diagram for small nuclear systems. 

The most troubling point in this otherwise elegant pic- 
ture is summarized by the question: where is the vapor? 
Does the nuclear system truly present itself at some time 
like a mixed phase system with the vapor bei ng s ome- 
how restrained (as in lattice gas calculations [_]]), ei- 
ther statically or dynamically in contact with the liquid 
phase? If so, how is the nuclear vapor restrained? If not, 
what is the meaning of vapor pressure, when clearly the 
system is freely decaying in vacuum against no pressure? 

In this paper we present a physical picture of fragment 
production from excited nuclei that will show: (a) how 
one can talk about coexistence without the vapor being 
present; (b) why a simple thermometric equation such 
as E = aT 2 is a more pertinent thermometer for inter- 
mediate mass fragment production than empirical ther- 
mometers such as isotope thermometers; and (c) why an 
equilibrium description, such as the FMD, is relevant to 
the free vacuum decay of a multifragmenting system. 

To begin, let us imagine a liquid in equilibrium with 
its saturated vapor. At equilibrium, any particle evapo- 
rated by the liquid is restored on the average by the vapor 
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bombarding it. In other words, the outward evapora- 
tion flux from the liquid to the vapor is exactly matched 
by the inward condensation flux. This is true for any 
kind of evaporated particle. The vapor acts like a mir- 
ror, reflecting back into the liquid the particles which it 
is evaporating. One obviously can probe the vapor by 
putting a detector in contact with it. However, since the 
outward and inward fluxes are identically the same, one 
might as well put the detector in contact with the liquid 
itself. At equilibrium, the two measured fluxes must be 
the same. Therefore, we do not need the vapor to be 
physically present in order to characterize it completely. 
We can just as well study the evaporation of the liquid 
and dispense with our imaginary surrounding saturated 
vapor. The vapor need not be there at all. One speaks 
in these situations of a "virtual vapor", realizing that 
first order phase transitions depend exclusively upon the 
intrinsic properties of the two phases, and not on their 
interaction. Of course, if the vapor is not there to re- 
store the emitting system with its back flux, evaporation 
will proceed, leading to a cooling off of the system. This 
naturally suggests the study of the emission of fragments 
from an excited nucleus since the nucleus is a small drop 
of nuclear liquid evaporating in vacuum; i.e. emitting 
fragments in vacuum. 

We now specifically address fragment emission from an 
excited nucleus with a time-honored assumption which 
we do not justify other than through the clarification 
it brings to the experimental picture: we assume (just 
as in compound nuclear decay) that, after prompt emis- 
sion in the initial phase of the collision, the resulting sys- 
tem relaxes in shape and density and thermalizes on a 
time scale shorter than its thermal decay. At this point 
the excited nucleus emits particles in vacuum, accord- 
ing to standard statistical decay rate theory. In this pic- 
ture there is no surrounding vapor, no confining box, and 
there no need for either. By studying the outward flux 
of the first fragments emitted, we can study the nature 
of the vapor even when it is absent (the virtual vapor) 
because of the equivalence of the evaporation and conden- 
sation fluxes of a liquid in equilibrium with its saturated 
vapor. 

Quantitatively, the concentration ua{T) of any species 
A in the vapor is related to the corresponding decay rate 
Ra(T) (or to the decay width Ta) from the nucleus by 
matching the evaporation and condensation fluxes 

Ra(T) = - M T ) (v A (T)4o- mv (v A )) , (4) 

where va(T) is the velocity of the species A (of order 
(T/A) 1 / 2 ) crossing the nuclear interface represented by 
the cross section a- lnv (of order A^ 3 where A is the 
mass number of the evaporating nucleus). The temper- 
ature To of the equilibrated, excited nucleus when the 
first fragment is emitted can be estimated by the thermo- 
metric equation of a Fermi gas and the calorimetrically 
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FIG. 1: The effective temperature of a Fermi system with 
three exit channels (a, b, and n) is plotted as function of 
initial excitation energy for two cases: one where barriers 
B a and £?{, are large (crosses) compared to B n =6 MeV, and 
another where B a or Bb is similar (solid circles) to B n . The 
initial temperature as a function of initial excitation energy 
is shown by the open circles. 

mesaured excitation energy E* 

T = y/Wfa (5) 

allowing for a weak dependence of a on T [l3l . IhH , and 
remembering that the system is most likely still in the 
Fermi strong degeneracy regime (where the temperature 
is much less than the Fermi energy: T <ti £f)- 

This is the fundamental and simple connection be- 
tween Eq. @, the (compound nucleus) decay rate, and 
Eq. (I), the FDM. In the latter, one immediately rec- 
ognizes in the exponential the canonical expansion of 
the standard compound nucleus decay rate, namely the 
Boltzmann factor exp(—B/T) where B is the emission 
barrier which in Eq. JJ| is written with its surface factor 
isolated from all other components, e.g. Coulomb [Io| . 
symmetry, finite size |9(, etc. Thus, the vapor phase in 
equilibrium can be completely characterized in terms of 
the decay rate. 

The physical picture described above is valid instanta- 
neously, but not globally. The result of a global evapora- 
tion in vacuum leads to abundances of various species of 
emitted fragments that arise from a continuum of systems 
at different temperatures. This leads to complications in 
various thermometers: kinetic energy, isotope ratios, etc. 
One way to avoid this complication is to consider only 
fragments that are emitted very rarely so that, if they 
are not emitted first, they are effectively not emitted at 
all. In other words, we consider only fragments that by 
virtue of their large surface energy, have a high emission 
barrier. 

As an example of the effects of this complication on the 
isotope thermometer, consider a decaying system with 
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only three available exit channels a, b, and n with barriers 
B a , B b , and B n . For B n <C B a and B b we know that the 
probability of emission of particles of type b at a fixed 
temperature is approximately 



p b «exp[- [B b -B n )/T]. 



(6) 



Since the nucleus cools as particles are emitted, the total 
emission probability of particles of type b from a nucleus 
at T is 

A similar expression holds for P a . The ratio of Pb/P a is 



A a r o To/A " e-V'xdx 



(8) 



where A& = B b — B n and A a = B a — B n . The ratio Pb/ P a 
can also be used to extract an effective temperature T e g 



P b /P a = exp [- (B b - B a ) /T eS ] 



(9) 



This sequential aspect is of course lost in equilibrium 
models which, at best, are corrected for the secondary 
decay of the primary fragments, but not for the source 
decay. 

A comparison of T e g and To is given in Fig. [2 for dif- 
ferent values of B b and B a . The case where B a and B b 
are large (crosses) gives effective temperatures very near 
to the initial temperature To (open circles) . When either 
B a or B b is near the barrier of the most probable channel 
(solid circles), the effective temperature is very different 
from the initial temperature. 

Thus, in order to justify the use of the initial Fermi 
temperatures one should choose exit channels with large 
barriers. This was done in the analyses leading to the nu- 
clear phase diagrams Q, , where fragments with charge 
Z < 5 were not considered. Under these conditions, the 
validity of Eq. Q is essentially guaranteed. The emis- 
sion rate can then be related to the vapor concentration 
and the phase diagram can be constructed. The tem- 
perature necessary for our purpose is To and not some 
average temperature determined from multiply emitted 
particles. 

The correctness of the thermometric relation given by 
Eq. JSJ can be tested "a posteriori" by verifying the lin- 
earity of the Arrhenius plots In tla vs. 1/T |l5|and of the 
socalled Fisher plots of ua scaled via Eq. Q 1, 2J. This 
linearity, extending over many orders of magnitude for a 
variety of fragments, is the strongest test yet of a Fermi 
gas thermometric relationship. In fact one can reverse 
the problem and determine the thermometric relation- 
ship up to rather high excitation energies by the require- 
ment that it leads to a linear Arrhenius and Fisher plots. 
This is an important, novel insight into the problem of 
temperature in excited nuclei. 
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FIG. 2: The mean emission times (in fm/c) of fragments 
with atomic number 4 < Z < 9 are plotted (solid symbols) 
versus inverse temperature for the reaction 7r+Au at 8 GeV/c 
[i6lll7|l . The average yields of the same fragments are plotted 
versus 1/T (solid symbols). The line represents a Boltzmann 
fit to the fragment yields. This same line has been superim- 
posed (shifted) on to the emission times. 



Further evidence for our physical picture of an excited 
nucleus evaporating fragments can be found in the mean 
emission time of fragments. To see this, we first construct 
an Arrhenius plot (tia{T) versus 1/T) of the abundances 
of observed fragments as a function temperature such 
as the one in Fig. [21 which shows the reciprocal of the 
average abundances 1 / (n) for fragments with 4 < Z < 9 
measured by the ISiS collaboration in the reaction 7r+Au 
at 8 GeV/c [IE E3l 

The slope of the Arrhenius plot is the effective "bar- 
rier" Ba for the emission of the fragment. This is evident 
in Eq. Q which gives 



n A (T) cx T A (T) cx exp(-B A /T). 



(10) 



Equation @ also shows that the same barrier and the 
same Boltzmann factor determine the mean emission 
time t of a fragment since 



tA(T) 



Ta(T) 



cx exp(B A /T). 



(11) 



Such a time £a(T) is the reciprocal of Ta{T). 

Therefore, the same Arrhenius plot with the same bar- 
rier should describe both the temperature dependence of 
the abundances and of the times. This is exactly the case 
shown in Fig. El The ISiS collaboration has measured the 
yields (open symbols) H| and the mean emission times t 
(solid symbols) 0, Il7j of fragments with 4 < Z < 9 
as a function of excitation energy. These energies can 
be translated into a Fermi gas temperature 0] as dis- 
cussed above. A Boltzmann fit to the yields is shown 
by the solid line. That same line has been superimposed 
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(shifted) onto the emission time data which it describes 
very well. Interestingly, the two different observables and 
their energy dependence are described by equations ljlTJ|) 
and Ijllfl with the same barrier. 

All that has been written above holds exactly for low 
excitation energies, thus ordinary compound nuclear de- 
cay is directly relevant to the liquid-vapor phase transi- 
tion. We can then find further evidence for our physical 
picture by scaling known low energy compound nucleus 
fragment yields [18J] according to the FDM. 

In compound nucleus (and multifragmentation) exper- 
iments we measure the fragment yields Ya(T) (the num- 
ber of fragments in events with a given excitation energy 
or temperature divided by the total number of events at 
that excitation energy or temperature) rather than the 
fragment concentrations ha (T) (the number of fragments 
of a given mass per unit volume). However, following the 
logic set forth above, we can relate the measured yields to 
the equilibrium fragment concentrations as follows: since 
the fragments in question are large (Z > 7) they are emit- 
ted from the compound nucleus of Aq nucleons first or 
not at all. Then for first chance fragments the yield is 
given by: 



Ya{T) = 



AtR A {T) 



AtR n (T) + Atj: A ^ n R A (T) 



(12) 



where At is the time duration of the measurement and 
R n (T) is the rate of neutron decay rate. Since neutron 
emission is much more probable than heavy fragment 
emission, in other words R n {T) ^S> Y^Ajtn Ra(T), the 
yield of a given fragment is 



Ya(T) 



Ra{T) 

Rn{T) 



Recalling that 



R n {T) 



T 
Ti 



exp 



T 



(13) 



(14) 



where B n « 8 MeV is the neutron binding energy. It 
follows from equations l|13|) and (|14J) that 



Ya(T) 



^exp ) n A lT) (r.i(D-lf7i, lv (r 1 )) 

, HA* 



o_ A -r ( B n + Am - c eA° 



where go' consolidates the constants from va(T), 
4amv(Vk) and ua(T). The measured first chance frag- 
ment yields are thus proportional to the concentration of 
the virtual vapor. 

The yields of charged fragments from the reaction 
64 Ni+ 12 C were fit to Eq. HSJ). These data were taken 
at the Berkeley 88-inch cyclotron using Ni beams with en- 
ergies between 6 and 13 AMeV suggesting Aq = 76, ver- 
ified experimentally |l8| and excitation energies of 0.96 
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FIG. 3: Results for the FDM-scaled yield distribution for the 
64 Ni+ 12 C compound nucleus decay data. See text for details. 



TABLE I: Results for fitting parameters 



16.8 MeV 



18.63 MeV 



cq 

~xl L39 L53 

go' 40 ± 20 300 ± 100 

A/i 2.43 ± 0.05MeV 3.01 ± 0.05 MeV 
Tc 10.6 ± 0.6MeV 14 ± 1 MeV 



AMeV < E* < 1.82 AMeV and temperatures (calcu- 
lated as in references [H H Q) of 2.9 MeV < T < 4.2 
MeV. These excitation energies are small and the frag- 
ment emission barriers are large compared to those of 
neutron evaporation; therefore there is little doubt about 
the validity of the thermometric relation of Eq. j£| . The 
mass of a fragment A was estimated from the measured 
charge of a fragment Z using the EPAX parameterization 
[l9) . Because the excitation energies are low, the effects 
of the secondary decay of the fragments is minimal. 

54 points for fragments with charges of 7 < Z < 15 
from the 64 Ni+ 12 C reaction were fit to Eq. (O with 
three free parameters: qo' , A/i and T c . The critical ex- 
ponents t = 2.209 ±0.006 and a = 0.63946 ±0.0008 were 
set to their standard three dimensional Ising values . 
As a test of the systematic errors the surface energy co- 
efficient co was set to its text book value 16.8 MeV 
from the Weizacker-Bethe semiempirical formula and to 
18.63 as suggested by more recent efforts |2l|. 

Results are shown in Fig. and given in Table ID The 
meaning of the value of the normalization qo' is unclear, 
though it must depend on the constants or proportional- 
ity in va{T), Uinv{VA) and tla{T). The effective chemi- 
cal potential Afj, value is due to the effects of finite size 
0, the energetic contributions from Coulomb and 
other nuclear effects (e.g. isospin, angular momentum, 
etc.) and the Coulombic contributions to <7i nv . It is 
assumed that Afj, is approximately constant for the 76 
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FIG. 4: The solid line shows the ratio of coe (the surface free 
energy coefficient) for the two estimates of the nuclear binding 
energy used as a function of temperature. Dashed lines show 
the error arising from the errors on the T c values returned by 
the fitting procedure. Dotted vertical lines show the temper- 
ature range covered by the compound nucleus experiment. 

nucleon system over the small temperature range con- 
sidered. The value of the critical temperature between 
10.6 ± 0.6 and 14 ± 1 MeV (whose value is essentially 
dictated by the choice of the surface energy coefficient 
of nuclear matter) can be interpreted as the critical tem- 
perature of bulk, neutral, symmetric nuclear matter since 
the effective chemical potential Afj, absorbs all but sur- 
face effects. 

The systematic error of the fit parameters in this anal- 
ysis is sizable. This is not unexpected since the range 
in temperature covered by the compound nucleus exper- 
iment is small compared to either value of T c thus prvid- 
ing a limited range over which to test the dependence of 
the fragment yieds on cqE (the surface free energy coeffi- 
cient). The solid line in Fig.0]shows the ratio of cqE for 
the two estimates of the nuclear binding energy used as a 
function of temperature, the dotted vertical lines showing 
the temperature range covered by the compound nucleus 
experiment. The value of this ratio is nearly constant for 
the temperatures considered here. The difference of the 
ratio from unity is accounted for by the differing values 
of go' and A//. Without a larger lever arm in tempera- 
ture and a better understanding of go' an d this level 
of systematic error is unavoidable. 

In conclusion, we have shown that: a) first chance 
emission from a hot thermalized source is a natural ex- 
planation for the scaling of experimentally observed in- 
termediate mass fragment yields according to Eq JQl; 
b) The simple relationship between emission rates and 
saturated vapor concentration neatly eliminates the un- 



physical assumption of a physical vapor in equilibrium 
with a liquid; c) Fermi gas like temperatures associated 
with first chance emission find their natural explana- 
tion and are justified by scaling plots such as Fig. 03 d) 
The emission rate picture explains the identical slopes 
in the Arrhenius plots for the mean multiplicities and 
the mean emission times shown in Fig. [5] e) Finally, 
the fragment emission from an experimentally certified 
low energy compound nucleus reaction scales according 
to Eq. and leads among other things to a value of 
the infinite nuclear matter critical temperature between 
10.6 ± 0.6 and 14 ± 1 MeV in agreement with theoretical 
estimates jH HI |H El El ■ 

This work was supported by the US Department of 
Energy. 



[1] J. B. Elliott et al, Phys. Rev. Lett. 88, 042701 (2002). 
[2] J. B. Elliott et al, Phys. Rev. C 67, 024609 (2003). 
[3] J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995). 
[4] A. Chernomoretz et al, Phys. Rev. C 64, 044605 (2001). 
[5] F. Gulminelli et al, Phys. Rev. Lett. 82, 1402 (1999). 
[6] M. D'Agostino et al, Phys. Lett. B 473, 219 (2000). 
[7] M. E. Fisher, Physics (N.Y.) 3, 255 (1967). 
[8] M. E. Fisher, Rep. Prog. Phys. 30, 615 (1967). 
[9] L. G. Moretto, K. A. Bugaev, J. B. Elliott, R. Ghetti, 

J. Helgesson, and L. Phair, Phys. Rev. Lett. 94, 202701 

(2005). 

[10] L. G. Moretto et al, Phys. Rev. C 68, 061602 (2003). 
[11] C. M. Mader et al, Phys, Rev, C 68, 064601 (2003). 
[12] J. B. Elliott, L. G. Moretto, and L. Phair, Phys. Rev. C 

71, 024607 (2005). 
[13] K. Hagel et al, Nucl. Phys. A 486, 429 (1988). 
[14] A. H. Raduta et al, Phys. Rev. C 55, 1344 (1997). 
[15] L. G. Moretto et al, Phys. Rep. 287, 249 (1997). 
[16] L. Beaulieu et al, Phys. Rev. Lett. 84, 5971 (2001). 
[17] L. Beaulieu et al, Phys. Rev. C 63, 031302 (2001). 
[18] T. S. Fan et al, Nuc. Phys. A 679, 121 (2000). 
[19] K. Siimmerer et al, Phys. Rev. C 42, 25462561 (1990). 
[20] K. S. Krane, Introductory Nuclear Physics (Wiley, New 

York, 1988). 

[21] W. D. Myers and W. J. Swiatecki, Phys. Rev. C 63, 
034318 (2001). 

[22] B. Friedmann and V. R. Pandharipande, Nucl. Phys. A 

361, 502 (1981). 
[23] H. R. Jaqaman et al., Phys. Rev. C 27, 2782 (1983). 
[24] H. R. Jaqaman, A. Z. Mekjian, L. Zamick, Phys. Rev. C 

29, 2067 (1984). 
[25] S. Levit and P. Bonche, Nucl. Phys. A 437, 426 (1985). 
[26] N. K. Glendenning, L. P. Csernai, J. I. Kapusta, Phys. 

Rev. C 33, 1299 (1986). 
[27] H Mueller and B. D. Serot, Phys. Rev. C 52, 2072 (1995). 
[28] J. N. De, B. K. Agrawal, S. K. Samaddar PRC 59, Rl 

(1999). 



